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ABSTRACT 


A simple strength of materials analysis was developed for a double- 
cantilever beam (DCB) specimen to account for geometric nonlinearity effects 
due to large deflections and T-tabs. A new DCB data analysis procedure was 
developed to Include the effects of these nonlinearities. The results of the 
analysis were evaluated by DCB tests performed for materials having a wide 
range of toughnesses. The materials used in the present study were T300/5208, 
IM7/8551-7, and AS4/PEEK. 

Based on the present analysis, for a typical deflection/crack length 
ratio of 0.3 (for AS4/PEEK) , T-tabs and large deflections cause a 15 percent 
and 3 percent error, respectively, in the computed Mode I strain energy 
release rate. Design guidelines for DCB specimen thickness and T-tab height 
were also developed in order to keep errors due to these nonlinearities within 
2 percent. 

Based on the test results, for both hinged and tabbed specimens, the 
effects of large deflection on the Mode I fracture toughness (^j^) were almost 
negligible (less than 1 percent) in the case of T300/5208 and IM7/8551-7, 
however, AS4/PEEK showed a 2 to 3 percent effect. The effects of T-tabs on 
were more significant for all the materials with T300/5208 showing a 5 
percent error, IM7/8551-7 a 15 percent error, and, AS4/PEEK a 20 percent 
error . 

Key Words : Double cantilever beam, delamination, fracture toughness, 

composite, large deflection, geometric nonlinearity, strain-energy release 
rate , loading tabs . 
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NOMENCLATURE 


a crack length including end correction 

effective crack length for DCB with large deflections 

effective crack length for DCB with T-tabs 

effective crack length for DCB with T-tabs and large deflections 

^exp’ ^i Pleasured crack length 

a end correction to account for DCB crack tip rotation 

o ^ 

b DCB specimen width 

C^ measured initial compliance for DCB specimen 

d tab height from ba¥e to loading point 

d f 6, measured displacement in a DCB test 

E^^ longitudinal composite modulus 

E ^2 transverse composite modulus 

E^ transverse modulus of resin-rich layer 

composite shear modulus 

strain-energy release rate 
e t 

Gj strain-energy release rate for DCB with T-tabs 

e 1 1 

Gj Gj for DCB with T-tabs and large deflections 

Gj^ Mode I fracture toughness 

h thickness of DCB arm 

H total tab height from loading point to center of DCB arm 

H factor used in calculation of G^ 

Ic 

I moment of inertia 

k foundation spring constant 

M bending moment at a section x from the crack tip 
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bending moment at crack tip 

resisting moment due to T-tab rotation 

applied load 

DCB specimen thickness 

thickness of resin-rich layer 

deflection of the beam at any distance x from the crack tip 
distance from the crack tip 
deflection of DCB arm 

crack length shortening due to large deflections 
crack length shortening due to T-tab rotation 
angle of rotation at DCB loading point 


INTRODUCTION 


The double cantilever beam (DCB) specimen is a popular specimen for 
determination of composite Mode-I interlaminar fracture toughness (G^^) . DCB 
specimens usually consist of several unidirectional plies layed-up with a thin 
insert at the raidplane (Fig. 1) to serve as a starter crack. Load is applied 
either through metal hinges [1 - 8] or metal T-tabs [9 - 16] bonded to the end 
of the specimen. The crack length is measured along the specimen edges. The 
measured crack length, applied load, and load-point deflection are then used 
to compute by a data analysis method which is usually based on linear 

beam theory [1,2,5-8,12-18], 

In general, the DCB specimen is designed to limit deflections to the 
geometrically linear range [1,2]. However, with the advent of tough resin 
composites, there is an increased possibility of encountering large 
deflections for DCB specimen thicknesses that have been conventionally used. 
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Large deflections cause an effective shortening of the crack length (Fig. 

2(a)) which leads to errors in the computed values if DCB data are 

analyzed using linear beam theory assumptions. In such Instances, for a 
chosen DCB specimen design, an estimate of the effect of large deflections on 
the computed values can be made by the analyses of references 3 and 4, 

Alternatively, the design criteria suggested in references 1 and 2 can be used 
to limit DCB deflections to the geometrically linear range. 

As mentioned earlier, the DCB specimen is loaded either through bonded 
hinges or T-tabs. The use of T-tabs shifts the line of action of the load due 
to the rotation of the DCB arras and leads to an effective crack- length 
shortening (Fig, 2(b)). The use of hinged tabs drastically reduces this 
effect. The effective crack- length shortening in the case of T-tabs increases 
with load [9-11] and leads to a geometrically nonlinear problem. The analyses 
in references 9 and 11 account for the geometric nonlinearities resulting from 
both T-tabs and large deflections. Wang et. al [9] and Williams [11] state 
that the effects of loading tabs significantly contribute to the strain-energy 
release rate in a DCB specimen. 

Hinged DCB specimens are, therefore, preferable to T- tabbed specimens. 
T-tabs may be required for some tough materials, however, because of the 
larger critical loads. Although the effect of end rotation on computed G^ 
values could be accurately accounted for by the analyses in references 9 and 
11, these analyses are quite complicated and tedious to use. Conversely, T- 
tabs could be designed to minimize the effects of end rotation. However, 
there are no design guidelines available for DCB specimen tabs. 

The purpose of this paper is first to present a simple strength of 
materials analysis to account for geometric nonlinearities resulting from T- 
tabs and large deflections. Next, DCB design guidelines are presented to 
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minimize the effects of these geometric nonlinearities. Then, a data analysis 
procedure is developed to compute Finally, DCB test results are 

presented to quantify these effects for commonly used composites over a wide 
range of toughnesses . 

MATHEMATICAL ANALYSIS OF DCB SPECIMEN WITH T-TABS 

The analysis of a DCB specimen with T-tabs is presented in the following 
sections. The DCB arm is first idealized as a beam on an elastic foundation. 
Equations for the load-deflection behavior are then derived for the case of 
small deflections followed by derivations for the case of large deflections. 
Next, expressions for the strain-energy release rate are derived for a DCB 
with T-tabs. Finally, design guidelines are developed to ensure that DCB 
deflections are in the geometrically linear range. 

The DCB specimen arm can be represented as a beam on an elastic 
foundation [18], The deflection of such a beam with an elastic foundation can 
be approximated by the deflection of a cantilever beam with an additional end 
correction length of a^ (see Appendix). Thus, for the DCB specimen shown in 
figure 3, each arm of the DCB can be idealized as a cantilever beam of length 
a given by 


a 


a 


exp 


+ a 

o 


( 1 ) 


where a is the crack length measured during testing (see Figs. 1, 3 and 
exp 

4) and a^ is given by (see Appendix): 
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= h V(E^^/6 E 22 ) 


( 2 ) 


a 


where h Is the thickness of each arm; and E ^2 the longitudinal 

and transverse moduli, respectively. 


Small Deflection Analysis 


Consider a DCB specimen that has bonded T-tabs of height d (Fig. 3). 
The application of load P causes an end rotation 0 which results in a 
resisting moment (Fig. 3) given by 


» P H sin(^) (3) 

where, H =- d + (h/2) (4) 

Note that the tabs are assumed to be rigid in the present analysis. Thus, a 
DCB specimen with load P applied through end tabs and a measured crack 
length of %xp analyzed as a cantilever beam of length a (Eq,(l)) 

with load P and a resisting moment (see Fig. 3) . Note that is a 

function of $ which in turn is a function of P leading to a geometrically 
nonlinear problem. 

The cantilever beam in figure 3 can be analyzed by first writing the 
moment at any point x along the beam and then using the moment -curvature 
relationship for a beam. The moment at any point x along the beam is, 

M(x) =- P (a - x) - P H sin($) (5) 


6 



The moment -curvature relationship for the beam is given by [19], 




( 6 ) 


where I is the moment of inertia of the beam cross-section and w is the 
beam deflection at any point along the beam. Note that the flexural rigidity 
of the tabs is neglected in the above equation. Integrating and applying the 
appropriate boundary conditions yields expressions for the slope and 
deflection at the point A (see Fig. 3). The slope is given as 

2 

tan(e) - ( 2 i )(l - 2 -^sin(g)) (7) 

and the deflection 8 at point A is given by, 

A 

Note that equations (7) and (8) agree with the slope and deflection of a 
linear beam (with load P) when H is set equal to zero. 

As shown in figure 2, the loading tabs cause an effective shortening (A^) 
of the beam length. This shortening can be easily determined, by inspection of 
figures 3 and 4, as (H sin(^)). The effective length of a beam with T-tabs 
can, thus, be approximated by where, 

Set “ a (1 - -^stnO?)) (9) 
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will be used later in the derivation of Mode I 


This effective length 

strain-energy release rate for a DCB with T-tabs, Note that a^^ from the 

above equation can be used in the deflection equation of a linear cantilever 

beam to approximate the deflection of a linear beam with tabs; however, that 

would give the deflection at point A' (see Fig, 3) on the beam and not at 

point A. Note that, A' is not a fixed point on the beam and does not 

provide a good reference for beam deflection and thus the point A is used as 

a reference for beam deflection. The deflection at point A given by 

equation (8) can also be related to the measured (load-point) displacement in 

a tabbed DCB test. This measured deflection d (see Fig. 4) can be 

exp ^ 

expressed as, 


"^exp " 2 ^ ” (cos(^) - D) (10) 

The load-deflection behavior of a DCB with tabs can now be examined by using 
equations (7) , (8) and (10) . 

Figure 5 shows a plot of DCB load verses load-point deflection. Both the 
quantities are presented in a nondiraensional form. The dash- dotted curve 
represents the load-deflection behavior of a linear DCB with tabs for a H/a 
ratio of 0.3 which represents a 25.4 mm tab and a crack length of 85 mm. Note 
that, the crack lengths used in a DCB test are usually between 50 mm and 120 
mm. For deflection/length ratios of less than 0,05, the DCB with tabs closely 
follows the simple linear DCB with hinges (dotted line). However, for larger 
de^p/a ratios, the DCB with tabs departs considerably from linear beam 
behavior. For a DCB specimen (3 mm thick) made from AS4/PEEK, initial crack 
extension occurs at a d /a ratio of about 0.3 under static loading [15]. 
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At this ratio there is a difference of 18 percent between a DCB with tabs and 
a DCB with hinges. 


Large Deflections Analysis 


The effects of large deflections on DCB response has been studied in 
references 3, 4, 9, and 11. DCB specimens with hinges were considered in 
references 3 and 4, while tabbed DCB's were analyzed in references 9 and 11, 
The analyses in references 9 and 11, however, are very complicated and do not 
separate the effects of T-tabs and large deflections. The present study uses 

a simple strength of materials approach, similar to that in reference 4, to 

analyze a DCB with tabs undergoing large deflections. 

As discussed earlier, large deflections cause an effective shortening of 
the crack length (see Fig. 2). This shortening can be derived by considering 
the x-axis strain at the midplane of the beam. The strain c at the 

X 

midplane is given by [20] 

'x - S + (1/2 x|^)2 (11) 

where u and w are displacements in the x and y directions, 
respectively. The midplane strain is assumed to be zero [4]. Thus, the 

variation of u with respect to x is given by 

d\x/dx * -l/2(5w/5x)^ (12) 


The term (9w/3x) can be obtained from equations (5) and (6) for a DCB with 
tabs and is given by 
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dx 


( P )(^2 . (3 . - 2 X H sln(6»)) 

2 I 


(13) 


Integrating equation (12) along the beam length after substituting for dw/dx 
from equation (13) gives the total shortening of the beam in the x- 

direction as 


J <g) dx 

0 


(14a) 


1 . e . 


— ^ - (1/3) ( )^((l/5) - (5/8)(H/a)sin«» + (1/2) (H/a)^sin^((?)) 

a 1 

(14b) 

The effective crack length ® beam with large deflections is thus 


el 


a (1 - (A^ /a)) 


(15) 


Using equations (9) and (15) it is clear that the effective crack length for a 
DCB arm with tabs undergoing large deflections can be written as 


elt 


a (1 - 


H 


■sin((?) - /a)) 


(16) 


This effective crack length for a DCB arm with T-tabs undergoing large 
deflections will be used later in the derivation of Mode I strain-energy 
release rate. Note that, if the effective crack length 1® used in the 

deflection equation of a linear beam in order to account for the effects of 
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tabs and large deflections, then that would give the deflection at point A' 
and not at point A, 

The deflection, at point A, of a DCB arm with tabs that is undergoing 
large deflections can be written, using equations (8) and (15) as 

3 

«A - ( 3"e^^ i )( 1 - 1-5 ^sin(«))(l - /a))^ (17) 

The nonlinear load-deflection response of a DCB with tabs can be plotted using 
equations (7), (10), (14), and (17) and is shown in Fig. 5 (solid curve). At 
a typical ratio of 0.3 [15] and a H/a ratio of 0.3, the nonlinear 

tabbed DCB (solid curve) is 6 percent above the linear tabbed DCB (dash-dotted 
curve) and 25 percent above the linear hinged DCB (dotted curve) . The load- 
deflection curve for the DCB with tabs and large deflections (solid curve) 
obtained by the present analysis agrees well with that obtained by the more 
accurate analysis of reference 9 (short dashed curve) for ratios of 

less than 0.4. However, for ^gj^pA ratios of 0,4 and greater there is more 
than 6 percent difference between the present analysis and that of reference 
9. The effect of large deflections in the absence of T-tabs can be examined 
by using equations (14) and (17) and substituting H — 0. This leads to the 
long dashed curve in figure 5 for a hinged nonlinear DCB. For a t 3 T>ical 
de^p/a ratio of 0.3, the effect of large deflections in a hinged DCB is only 
3 percent. 


11 



Strain-Energy Release Rate Analysis 


Based on the deflection equation for a linear DCB with tabs (equation 
(8)) and its similarity to the linear cantilever beam equation, it is possible 
to derive a simple equation for the strain- energy release rate. As discussed 
earlier, the geometric nonlinearities associated with large deflections and T- 
tabs lead to an effective shortening of the beam. Equations (9), (14), and 
(15) give a good estimate of this shortening as a function of the applied 
load. 

The strain-energy release rate for the DCB is given by [4] 


G 


I 


M 


b I 


( 18 ) 


where b is the width of the specimen and is the bending moment at the 

crack tip. For a linear hinged DCB, is given simply by (P a). Equation 

(18) is valid for a linear hinged DCB, however, it has been shown in 
references 4 and 11 that it can be used for a DCB with geometric 
nonlinearities if the moment is calculated by taking into account the 
shortening of the DCB arm. Thus if M is given by (P a^^) for a linear 
tabbed DCB (see equation (9)) and by (P for a tabbed DCB with large 

deflections (see equation (16)) then equation (18) can also be used for these 
cases. The strain-energy release rate for a linear DCB with tabs can, 
therefore, be written as 


12 


2 2 

°r- 

Note that a in the above equation is the sum of the measured crack length 
a^xp correction a^ (see equation (1)), The angle 0 is a 

function of the applied load (see equation (7)), 

For a DCB with tabs undergoing large deflections, the expression for 
is given by, 


2 2 

.elt ,_P a . 

’I ■ 


■^sin(«) - (A^ /a)) 2 


( 20 ) 


where /a is the shortening in the crack length due to large deflections 
and is given by equation (14). The effects of both large deflections and T- 
tabs on the strain- energy release rate in a DCB can now be studied for a range 
of d /a ratios by using equations (7), (8), (10), (14), (17), (19), and 
( 20 ). 

2 

Figure 6 shows the normalized strain-energy release rate (G^ba /E^^I) as 

a function of the normalized DCB load-point deflection ^ 

ratio of 0.3. For d /a ratios of less than 0.3, the short dashed curve 

exp^ 

for a linear DCB with tabs differs by less than 4 percent from the results, 

for a linear DCB with tabs, from reference 10 (dash-dotted curve). At d /a 

exp' 

ratios greater than 0.3, the agreement is not very good; but for such high 
ratios there will be large deflections in the DCB specimen; and, one should 
use the results for the large deflection case. The solid curve in figure 6 
corresponds to a DCB with tabs and large deflections (equation (20)) and 
compares very well with the result from reference 9. For a typical d^xp/^ 
ratio of 0.3 [15], the present solution differs by only 2 percent from the 
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more accurate solution of reference 9. Also, for ratios greater than 

0.3 there Is a very good correlation between the present solution and that of 
reference 9. 

Figure 7 shows the effect of tab height for a range of H/a ratios 

(based on equation (20)). The solid curve corresponds to a nonlinear hinged 

DCB and was obtained from equation (20) by substituting H/a equal to zero. 

For a typical ratio of 0,3, the results for the nonlinear hinged DCB 

differ by only 3 percent from those for the linear hinged DCB (dotted curve) . 

The results for the linear hinged DCB were obtained by substituting H/a 

equal to zero in equation (19). The difference between the linear hinged DCB 

and the nonlinear tabbed DCB increases with increasing H/a ratios. 

The percentage variation from the linear beam theory assumptions for DCB 

specimens, that have T-tab and/or large deflection effects, are svunmarized in 

figure 8 for a range of d^^p/a ratios and a H/a ratio of 0.3. For a 

nonlinear hinged DCB, the percentage error in using linear beam theory 

analysis is less than 5 percent for d /a ratios of less than 0.4 

exp 

However, for a DCB with tabs (H/a -=0.3) undergoing large deflections, there 
could be errors as high as 18 percent (for a d^^^/a of 0.3) if one uses 
linear beam theory assumptions. For a linear DCB with tabs (H/a - 0.3) there 
would be a 15 percent error, for ad /a of 0.3, if one did not account for 

6Xp 

the effects of the tabs. 

Guidelines for Minimizing Geometric Nonlinearity Effects 

The Gj expressions in the previous section can also be used to derive 
design criteria for DCB specimens with T-tabs and large deflections to ensure 
that DCB deflections are in the geometrically linear range. 
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DCB Specimens with Large Deflections 


Consider a hinged DCB specimen which is undergoing large deflections. 

The expression for in this case can be derived from equation (20) by 

substituting H equal to zero and the corresponding expression for from 

equation (14) 

<=1 - ( 21 ) 

Notice, that the first term in the above equation corresponds to the for 

a linear DCB and the second term corresponds to a correction attributable to 
large deflections. If one desires to limit large deflection effects to, say, 
less than 2 percent, then the following inequality should hold true: 

2 

(1 - (l/15)(r^)2)2 2, 0.98 (22) 

hi ^ 

Using the above inequality and equation (21) and substituting for I 
appropriately, an inequality for laminate thickness t can be written as 


t > 8.65 }(Gj^ aVii) <23) 

where is the Mode I fracture toughness for the material being tested and 
a is the crack length of the DCB specimen. The laminate thickness used for a 
DCB specimen should satisfy the condition above in order to ensure less than 2 
percent effects due to large deflections. The above inequality could also be 
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used to find the crack length a for which errors due large deflection will 
be less than 2 percent. If the value for crack extension from the 

insert is of Interest then, in order to keep large deflection effects below 2 

percent, the insert length should be less than (0.04 t for a 

given laminate thickness t. 

DCB Specimens with T-tabs 

The expression for , in the case of a DCB specimen with tabs, was 

derived earlier (see equation (19)). In equation (19), the first term 
corresponds to the G^ for a linear DCB and the second term corresponds to a 
correction attributable to the effects of T-tabs. If one desires to limit T- 
tab effects to less than 2 percent, then the following inequality should hold 
true : 


(1 - ^sin(i9))^ > 0.98 

a 


(24) 


Using the above inequality and equations (7) and (19) and substituting for I 
appropriately, an inequality for total tab height H can be written as 


H < 0.01 7(0.0434 + a^) (25) 

Note that H is the total tab height and is given by equation (4) . The tab 
height used in a DCB specimen with tabs should satisfy the condition above in 
order to ensure less than 2 percent effects due to the tabs. For a hinged 
specimen H corresponds to the height of the hinge axis of rotation above the 
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centerline of the DCB arm. For DCB specimens that are made thicker in order 
to minimize large deflection effects, the distance H even for hinge tabs 
could be significant and equation (25) should be used to check for this tab- 
like effect. 


DCB TESTING AND DATA ANALYSIS 


The results of the analysis in earlier sections suggests that the effects 
of geometric nonlinearity associated with T-tabs, on the computed Gj. values, 
can be as high as 18 percent, while the effects of large deflection will 
usually be less than 5 percent. In order to illustrate these results, 
specimens made from unidirectional composites having a wide range of Mode I 
fracture toughness values were tested both as hinged DCB's and as tabbed 
DCB's. The DCB data was analyzed using three different procedures. First, a 
data analysis procedure based on the present analysis was used to account for 
large deflection and T-tab effects. Next, the well known Berry procedure [17] 
based on linear beam theory was used. Finally, the area method was used to 
provide an accurate reference for comparing average G^^ values. 

Materials and Specimens 

Specimens were about 3 mm thick and were cut from unidirectional, 24 ply, 
panels with a Kapton film (0.0127 mm thick) crack starter at the midplane. 
Three different panels were made from T300/5208, IM7/8551-7 [21], and AS4/PEEK 
prepreg according to the manufacturer's instructions. After curing, the 
panels were cut into 152 mm by 25 mm DCB specimens. Only two specimens of 
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each material were tested since the focus of the present tests was to 
illustrate the analysis and compare different data analysis procedures. For 
each material, one of the specimens was bonded with aluminum alloy hinges and 
another was bonded with aluminum T-tabs that were 25 mm in height (d - 25 nun, 
see Fig. 3). Edges of the specimens were painted with white water-based 
typewriter correction fluid and marked at increments of 2.5 mm for visual 
crack measurements. 


Test Procedure 

Static tests were performed under displacement control in a screw-driven 
machine at a constant cross-head rate of 0.0085 mm/s. Load-deflection data 
was collected through a digital data acquisition system and was also plotted 
on a x-y plotter. Travelling microscopes were used to monitor crack length 
along both edges of the specimen. Load and deflection were noted for the 
initial crack extension from the insert and then at every 2.5 mm of crack 
extension indicated by the marks on the specimen edges. After the crack had 
extended about 12.7 mm the specimen was unloaded. The specimen was reloaded 
and the crack was extended another 12.7 mm while the load, deflection, and 
crack length were monitored. A third loading-unloading cycle was conducted in 
a similar manner for another 12.7 mm of crack extension. Typical load- 
deflection plots are shown schematically in figure 9, The geometric 
nonlinearities associated with T-tabs and large deflections cause an upwardly 
concave curve with a raonotonically decreasing compliance. The initial 
compliance, denoted by , will be used later in the data analysis 

procedure . 
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Data Analysis Procedures 


As mentioned earlier, data analysis was performed using three different 
procedures. This section describes the present analysis to account for the 
effects of large deflections and T-tabs, Berry's [17] procedure and the area 
method. 

Present Data Analysis Procedure 

The present procedure uses the initial compliance (Fig. 9), the load, 

deflection, and crack length data, and the expression given by equation 

(20) to determine There are, however, two unknowns in equation (20) 

that need to be determined before that equation can be used to determine , 

Ic 

These two unknowns are the flexural stiffness the angle $ at the 

end of the DCB arms . 

The quantity Eiil is directly evaluated by compliance calibration in 
the present data analysis procedure. However, as Indicated in figure 9, the 
compliance of a DCB with geometric nonlinearities changes with load. Based on 
the analysis, the load-deflection plots In figure (5) also indicate a 
compliance that changes with load but for very small deflections (5/a < 0.05), 
the curves for the geometrically nonlinear cases coincide with the linear beam 
curve. Thus the effects of large deflections and T-tabs are minimal for the 
initial part of the load-deflection curve and the initial compliance can 

be approximated by linear beam theory assumptions as 
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c. 

1 


(26) 


■3 Ej^r 


3 

a 


The quantity in brackets can be determined from a least - squares fit on a log- 
log plot of versus a. This method was used to determine the 

present data analysis procedure. 

The above procedure eliminates the errors involved in using the values of 

and I in the computation of quantity ^ in equation (20) 

is actually an effective modulus which, in general, is not the same as the 

inplane modulus [1,2,13]. Furthermore, the in a composite could vary 

along the laminate thickness due to inhomogeneity caused by the manufacturing 

3 3 

process [22]. Also, the moment of inertia I (*= (1/12) b h ) contains h 
which could lead to large errors in the computed I due to even small errors 
in measuring laminate thickness. 

The angle d at the end of the DCB arms was determined by using the 
load-deflection relation derived earlier (equation (17)) together with the 
result from equation (10). The measured deflection can thus be written as 

3 

‘^exp “ - 1.5 -J-sin(^)] + 2 H (cos(^) - 1) (27) 


A shear correction term [8] given by (2.4 P a/(b h G^^)) could be added to 

the above equation where the shear modulus and is determined 

as discussed earlier. The measured load P, deflection d and crack 

exp 

length a (given by equation (1)) are substituted into equation (27) to 
determine the angle $ iteratively by a numerical scheme. The secant method 
was used in the present study. 
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Once the quantities and 9 are known, the for the material 

can be computed by 

2 2 

*=10 - - (Aj /a))2 (28) 

where (A^^/a) is given by equation (14). A shear correction term [ 8 ] given by 
2 2 

(1.2 P /(b h ^^^ 3 )) could be added to the above equation to account for shear 

deformation. Note that the term ( (H/a) sin(^ ) ) corresponds to a correction 

due to T-tabs and should be used only for a DCB with tabs. Also, the term 

(A^ /a) corresponds to a correction due to large deflections. 

Alternatively, the could also be calculated by first determining 

Efil , as described earlier, and then using figure 7 along with the measured 

d^xp/^ ratio at crack extension and the H/a ratio to find the corresponding 
2 

value for (G^b a Knowing b, a, and one can find Gj and then 

calculate G^^ by adding the shear correction of equation (28) to G ^ . 

Berry's Method 

The mode I fracture toughness can be determined by Berry's method [17] 
which uses the load, load-point deflection, and crack length data at the time 
of crack extension in the following equation [13]: 


where , 


n H 

"Ic “ 2 b 


N 

H > (1/N) S (-y^) 

i=L '*1 


(29) 


(30) 
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where and are the measured load, deflection, and crack length, 

respectively, and N is the total number of data points while n is the 
slope of the least- squares fit to the compliance-crack length data plotted on 
a log-log scale. This method has the advantage of not requiring as the 

input parameter. Also, since the load, deflection, and crack length are all 
used in the computation, it might account, at least in part, for the effects 
of geometric nonlinearity. However, Berry's method is based on linear beam 
theory and does not explicitly account for the effects of T-tabs and large 
deflections . 

The Area Method 


The load-deflection plots shown in figure 9 can be used directly to 
compute by accurately measuring the area enclosed by the loading- 

unloading curve and dividing it by the incremental, area created during crack 
extension. This method implicitly accounts for any geometric nonlinearities 
since it uses the actual load-displacement curves. However, it can only give 
average values and cannot be used to compute the G^^ at the onset of 

crack growth from the insert. 


RESULTS AND DISCUSSION 

The DCB data for the three different materials was used to compute G 

Ic 

for both the hinged and the tabbed specimens. The present data analysis 
method was used by including all the geometric nonlinearity terms and then 
also used by neglecting the effects of large deflections and T-tabs. Table 1 
shows a comparison of the G^^ values computed using the present method, 
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Berry's method and the area method. The values shown in the second through 
fifth columns are for the onset of dolamination from the insert while the 
values shown in the last two columns are average values for the first 12 . 7 mm 
of crack extension. The values calculated using Berry's method (BM) should be 
compared with those from the present method obtained by neglecting the effects 
due to tabs and large deflections (LB). In general, there is good agreement 
between these two methods and the slight differences can be attributed to the 
differences in the two methods of data analysis. The average values 

shown in the last two columns are subject to the effects of fiber-bridging and 
are therefore higher than those computed for the onset of delamination. 
However, a comparison of these two columns helps validate the results of the 
present analysis technique. 

Table 1 also shows values computed using the present data analysis 

method after separately accounting for large deflection (LDT) and T-tab (TB) 
effects. For both hinged and Cabbed specimens the effects of large deflection 
were less than 0.8 percent in the case of T300/5208 and IM7/8551-7, and about 
3 percent for AS4/PEEK. The effects of T-tabs were more significant for all 
the materials with T300/5208 showing a 5 percent effect, IM7/8551-7 a 15 
percent effect, and, AS4/PEEK a 20 percent effect. 

The effects due to large deflections and T-tabs could have been avoided 

for the three materials if the DCB thickness and T-tab height were selected 

according to the guidelines given by equations (23) and (25), For T300/5208, 

2 

the two equations yield (using = 37.5 J/m , =» 181.3 GPa, a - 64 mm) 

t > 1.08 mm and H < 16 mm. In the present study a laminate thickness of 3 
mm and H - 25.75 mm were used, thus, leading to negligible large deflection 
effects and small T-tab effects. For AS4/PEEK, the guidelines of equations 
(23) and (25) yield (using = 1622 J/m^, - 136.5 GPa, a - 64 mm) t > 
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3.14 mm and H < 4 mm. The laminate thickness used in this case was 3.23 mm, 
thus, explaining large deflection errors of about 2 percent (see Table 1). 

The tab height H used was 25.8 mm which is much higher than the recommended 
4 mm leading to the 20 percent errors due to T-tab effects. 


CONCLUDING REMARKS 


A simple strength of materials analysis was developed for a DCB specimen 
to account for geometric nonlinearity effects due to large deflections and 
tabs. A new DCB data analysis procedure was developed to include the effects 
of these nonlinearities. The results of the analysis were validated by DCB 
tests performed for materials having a wide range of toughnesses. The 
materials used in the present study were T300/5208, IM7/8551-7, and AS4/PEEK. 

The results of the present simple analysis compared very well with 
previously developed, more complicated analyses. Based on the analysis, for a 
typical deflection/crack length ratio of 0.3 (for AS4/PEEK) and a H/a ratio 
of 0.3, there could be a 19 percent effect due to T-tabs and a 6 percent 
effect due to large deflections on the DCB load-deflection response. The 
computed strain-energy release rates can be in error by 15 percent due to tabs 
and by 3 percent due to large deflections for the same deflection/crack length 
ratio of 0.3. In order to keep errors due to these nonlinearities within 2 

3 2 

percent the DCB specimen thickness should be greater than 8.65 ^ ^^11^ 

3 2 

and the total tab height should be les.s than 0.01 7(0.0434 t + a ). 

Based on the test results for both hinged and tabbed specimens, the 
effects of large deflection were almost negligible (less than 1 percent) in 
the case of T300'/5208 and IM7/8551-7, however, AS4/PEEK showed a 2 to 3 
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percent effect. The effects of T-tabs were more significant for all the 
materials with T300/5208 showing a 5 percent effect, IM7/8551-7 a 15 percent 
effect, and, AS4/PEEK a 20 percent effect. The average values computed 

using the present analysis compared well with those calculated using the area 
integration method. 

APPENDIX 

The end deflection of an orthotropic beam on an elastic foundation can be 
written in a form similar to that given for an isotropic beam in reference 
[18] and is given as 

6 - ( , I + 3/(Aa) -I- 3/(Aa)^ + 1.5/(Aa)^) (Al) 

3 E^^ I 

where, A - *J((0.25 k)/(Ej^^ I)) (A2) 

and k is the spring constant of the elastic foundation. In the case of a 
DCB specimen the elastic foundation consists of a thin resin rich layer (of 
thickness t^) that forms in between two plies and an orthotropic laminate 
layer (of thickness h) [23]. The combined stiffness of the two layers 
represents the foundation spring constant k. By assuming a constant strain 
in the resin rich layer and a linearly varying strain distribution in the 
laminate layer [23] the tranverse stiffnesses are given by (E^ b/t^) and 
(2 ^*22 b/h) , respectively (where b is DCB width) . Since the resin rich layer 
and the laminate layer are in series , the foundation spring constant is given 

by 
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k 


2 E22 h 


(A 3 ) 


For a very thin resin rich layer equation (A 3 ) can be simplified as 


k =« 2 E22 b/h 


(A 4 ) 


and from equation (A 2 ) 


1 /A = h V(E^3^/6 E22) (A 5 ) 

Now, equation (Al) can be approximated by replacing 1.5 in the last terra with 
1.0 as 


5 


<3 I><‘ * 


(A6) 


where 


% “ 1 /A (A 7 ) 

Note that ecjuation (A6) represents the load-point deflection of a cantilever 
beam of length (a + a^) and thus a beam on elastic foundation could be 
approximated by a cantilever beam with an "end correction" a^ which is given 
from equations (A 5 ) and (A 7 ) by 


h ^y(E^j^/6 ^22^ (A®) 
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Since the present analysis is primarily for static DCB tests, only the 
extensional stiffness of the foundation was considered while it's rotational 
stiffness was neglected. It has been shown in reference 24 that for static 
cases the rotational stiffness can be neglected but it should be included in 
analyzing dynamic cases . 
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Table 1:- Comparison of DCB test results. 


Material 



^ initial 


(G^ ) (J/ni' 

Ic average ' 

Hinged 


LB^ 


LDT^ 

Area 

ldt'^ 

T300/5208 

94 

107 

- 

107 

112 

119 

IM7/8551-7 

605 

629 

- 

634 

756 

670 

AS4/PEEK 

1564 

1593 

- 

1622 

1705 

1784 

T- Tabbed 







T300/5208 

96 

102 

107 

107 

85 

100 

IM7/8551-7 

519 

536 

618 

621 

858 

771 

AS4/PEEK 

1429 

1401 

1680 

1734 

1873 

1805 


^ BM - Berry's method [17], 

^ LB - Present method neglecting tab and large deflection effects, 
c 

TB - Present method considering only tab effects. 

^ LDT - Present method accounting for tab and large deflection effects. 


31 




32 


Figure 1 , - DCB 



(b) Effect of large deflections and T-tabs 




Double cantilever beam specimen 
with end tabs 



Figure 3.- Cantilever beam idealization for a DCB specimen with T-tabs. 





Figure 4.- Nomenclature for a DCB specimen with T-tabs. 


Nonlinear DCB (tabbed) 
Nonlinear DCB [Ref, 9] 



Figure 5.- Normalized load- deflection curves for hinged and Cabbed DCB 
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Figure 8.- Percentage variation in G for DGB specimens with geometric 





a. 
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Figure 9.- Typical load-deflection curves for DCB specimens with geometric 
nonlinearities , 
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